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ABSTRACT 

We introduce a novel technique for empirically understanding galaxy evolution. We 
use empirically determined stellar evolution models to predict the past evolution of 
the Sloan Digital Sky Survey (SDSS-II) Luminous Red Galaxy (LRG) sample without 
any a-priori assumption about galaxy evolution. By carefully contrasting the evolu- 
tion of the predicted and observed number and luminosity densities we test the passive 
evolution scenario for galaxies of different luminosity, and determine minimum merger 
rates. We find that the LRG population is not purely coeval, with some of galaxies 
targeted at z < 0.23 and at z > 0.34 showing different dynamical growth than galaxies 
targeted throughout the sample. Our results show that the LRG population is dynami- 
cally growing, and that this growth must be dominated by the faint end. For the most 
luminous galaxies, we find lower minimum merger rates than required by previous 
studies that assume passive stellar evolution, suggesting that some of the dynamical 
evolution measured previously was actually due to galaxies with non-passive stellar 
evolution being incorrectly modelled. Our methodology can be used to identify and 
match coeval populations of galaxies across cosmic times, over one or more surveys. 
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1 INTRODUCTION 

The most recent observational evidence points towards a 
model of structure and galaxy formation that is hierarchi- 
cal in nature: small fluctuations in the matter density field 
grow via gravity, with their dynamics being governed in 
detail by dark matter and dark energy. Baryons trace the 
dark matter and, in regions of sufficient gravitation al depth, 
they accumulate and form stars and galaxies (Whit e~fc ReesI 
1 19781 ). As the matter field continues to evolve dynamically - 
largely oblivious to this process - the stellar content of galax- 
ies grows via a combination of two modes: forming new stars 
from cold gas, and merging with other galaxies. 

LRGs are extensively used as cos mological probes (e.g. 
iReid et al.l l2010l : iPercival et afl |2010| ). and are also inter- 
esting from a galaxy evolution perspective, as they domi- 
nate the galaxy mass function at the massive end. For these 
reasons, the steUar and dynamical evolution of LRGs and 
early- type galaxies (ETGs) has been studied extensively (see 
iToieiro fc PercivafcOldlToieiro et~alll201ll - henceforth TU 
- for a summary and list of references). Studies that address 
one of these growing modes, however, traditionally assume 
a model for the other. In this paper we propose and apply 
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a new methodology that solves for these two modes inde- 
pendently. In brief, we use the fossil record of local galax- 
ies to predict number and luminosity densities at past red- 
shifts. The differences between the predicted and the ob- 
served quantities, under some simple assumptions, can be 
interpreted as a merger history. 

The idea of inferring the properties of galaxies at a 
different cosmic time from the one they are observed at, 
and comparing them to the in-situ properties of galaxies at 
those cosmic times has been proposed before. Most notably 
iDrorv k, Alvare3 (|2008l ) use the galaxy stellar mass function 
(GSMF) of galaxies in the FORS Deep Field d Prorv et all 
12005) and estimates of the instantaneous star-formation rate 
as a function of observed stellar mass to separate the evolu- 
tion of the GSMF in terms of its merging and star formation 
components. Our approach differs from theirs in terms of 
the observables (we focus on number and luminosity densi- 
ties, as opposed to the GSMF), but mainly in how we effec- 
tively link the gala:xies at different cosmic times. We rely on 
the fossil record of local gala:xies to reconstr uct their past 
stellar build-up, which we get from VESPA ([Toieiro et al.l 
2007) in a non-parametric way. iDrorv fc Alv arez (2008') use 
parametric functions for the evolution of the GSMF and for 
the instantaneous star-formation rate, which they integrate 
over past cosmic time to predict a mass build up due to 
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star-formation only. In other words, we use present-day in- 
formation to predict the past Universe, whilst they use past 
information to predict the present-day Universe. Our recon- 
struction of stellar mass build-up is less parametric, but ar- 
guably more dependent on the underlying stellar population 
synthesis (SPS) models. The two results are not comparable 
due to the sample of galaxies studied, but the approaches 
are similar in terms of philosophy. An a pproach more simi- 
lar to ours was very recently presented in lEskew fc Zaritskvl 
|20l3), where the authors use published estimates of the 
past star formation histories of three local group galaxies 
(obtained using resolved stellar populations) to infer where 
these galaxies would sit in popular diagnostic plots at previ- 
ous times in their cosmic history. The work we present this 
paper is similar in terms of concept, but is vastly different in 
terms of the size and type of galaxies used, as well as overall 
goal. 

This paper is organised as follows: in Section [2] we in- 
troduce the data set we use, in Section [3] we introduce our 
methodology and in Section|4]we define our observables. We 
present our results in Section O which we interpret and dis- 
cuss in Section [6l We summarise and conclude in Section [T] 
Throughout this pape r we use assume a WMAP7 cosmology 
ijKomatsu et al.ll201ll '). 



2 DATA 

The SDSS is a photometric and spectroscopic survey, un- 
dertaken using a dedicated 2.5m telescope in Apache Point, 
New Mexico. For details on t he ha rdwa re, software and 
data-reduction see lYork et all l|2000l ) and IStought on et al.l 
In summa ry, the survey wa s carried out on a mo- 
saic CCD camera (|Gunn et al.lll99^ ). two 3-arcsec fibre-fed 
spectrographs, and an auxiliary 0.5m telescope for photo- 
metric calibration. Photometry was taken in five bands: 
u, g, r, i and z, and Luminous Red Galaxies (LRGs) were 
selected for spectrosco pic follow-up according to the target 
algorithm described in lEisenstein" et al.1 (|200ll ). In this pa- 
per we analyse the lat est SDSS LRG sample (data release 7, 
lAbazaiian et al.ll2009l ). which includes around 180,000 ob- 
jects with a spectroscopic footprint of around 8000 sq. de- 
grees and a redshift range 0.15 < 2: < 0.5. 

The target selection was designed to follow a passive 
stellar population in colour and apparent magnitude space, 
and it targets LRGs below and above z < 0.4 with two 
distinct cuts. Cut I targets low redshift LRGs by using the 
following cuts: 



rp < 13.1 -HC||/0.3 
rp < 19.2 

cx < 0.2 (1) 
^r.p < 24.2 mag arcsec^ 

'psf f model > 0.3 

where the two colours, cn and c± are defined as 



rosian magnitudes for the apparent magnitude and surface 
brightness cuts. Cut II targets LRGs at z > 0.4 following: 



rp < 19.5 
cx > 0.45 - (g 



r)/6 

(ff-r) > 1.3-f 0.25(r -i) 
/ir.p < 24.2 mag arcsec^ 

'^psf r-inodel ^ 0.5 



(3) 



Two separate cuts are necessary as the passive stellar 
population turns sharply m & g — r ys r — i colour plane, 
when the 4000A break moves through the filters. This bend 
in the colour selection results in a broader colour range of 
targets around it. 



3 METHOD 

Our approach consists of tapping into the fossil record of 
low-redshift [z < 0.25) LRGs to measure their stellar evo- 
lution in terms of their past star-formation and metallicity 
history out to 2: = 0.45 (Til). As discussed in Til, using 
the fossil record allows us to explore past star-formation his- 
tories in a way that is completely decoupled from the sur- 
vey's selection function. Although the SDSS selection was 
designed to follow a passively evolving stellar population, it 
is not true that all galaxies that pass the selection criteria 
follow this model. 

With the full knowledge of how the selection function 
of the survey changes with redshift, we predict the number 
and luminosity densities of LRGs at higher redshifts based 
on low-redshift samples. In the absence of mergers, the pre- 
dicted and observed number and luminosity densities should 
match at high redshift. In Section [6] we show how, by mak- 
ing a few simple assumptions, the difference between the two 
quantities can be interpreted as a minimum merger history. 



3.1 Evolving galaxies back in time 

For each galaxy we consider the past history, and at which 
epochs it would have been included in the survey. To do this 
we need to consider the following: 

(i) the stellar evolution of the galaxy; 

(ii) the selection function of the survey; 

(iii) changes in the photometric errors with redshift. 



(2) 



C|| = 0.7(ff - r) + 1.2[(r - i) - 0.18] 
c± = {r~i)-{g- r)/4 - 0.18. 

Model magnitudes are used for the colours, and pet- 



For (i) we u s e the stellar evolution models for the LRGs 
of iToieiro et al.l (|201ll 'l. and refer the reader to that paper 
for full details of how these models were obtained. For com- 
pleteness, we present a summary in Section [3.21 

To account for |(ii)[ we run the evolution of each indi- 
vidual galaxy thro ugh the selection cuts of the LRG sample 
(see Section and [Eisenstein et al.|[2"00lh . between the red- 
shift of the galaxy and z — 0.45, and we construct a vector 
that tells us exactly when any individual galaxy would have 
been observed within the survey, were it to exist on our past 
light cone at all epochs. We call this vector Vobs{z) and we 
set it to unity if the galaxy would have been observed at 
that redshift, and to zero otherwise. 

To incorporate the changing photometric errors ( |(iii)[ ) 
into our colour modelling we take the following steps. First 
we compute the average photometric errors in colour as a 
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function of apparent r— band cmodel magnitude, using all 
LRGs in the sample. Second, for each galaxy we compute 
the g — r and r — i offsets with respect to their assigned 
model at the redshift of the galaxy - note that these are 
not expected to match as each cell used to compute the 124 
stacks has a finite width in colour and redshift. We then add 
a term to their predicted colour evolution which keeps the 
ratio between this offset and the typical observational error 
at that redshift constant throughout its whole evolution. For 
example, to predict the g—r colour of a galaxy observed with 
a redshift Zobs at any other redshift z > Zobs we compute 



[g-r]pred{z) = [g - r]^odei{z) + A[g - r] 



[r{z)] 



ag-r[r{Zabs)] 



(4) 



where [g — r\modei is given by the stellar evolution models, 
^ig-r] = [g-r]obs - [g - r]modi,i{Zobs) , and r{z) is the pre- 
dicted apparent r— band cmodel magnitude at that redshift. 
When z — Zobs this simply returns the observed colour, and 
the effect at z > Zobs is to correctly match wider areas of 
colour space at larger redshift with narrower cells at low red- 
shift. Note that if there was no change in photometric errors 
with redshift there would be no need to widen the cells. 



3.2 The stellar evolution models, K+e corrections 
and absolute magnitudes 

We follow iToieiro &: Percivall (|2010l ) to compute K-l-e cor- 
rections and obtain K-f e corrected rest-frame absolute mag- 
nitudes - see Section 2.1 of that paper for details - with 
two difference s. First, we do not ass ume the passive evolu- 
tion model of iMaraston et aD (|2009l ) but rather the stellar 
evolution models of Til, which computed directly from the 
fossil record of LRGs. This allows us to proceed without 
making any assumptions about this evolution - thus mak- 
ing the results scientifically very robust. Secondly, a conse- 
quence of not wishing to make unsupported assumptions to 
construct evolutionary models is that we cannot predict how 
any galaxy would look like at z < Zobs- We therefore K-l-e 
correct all galaxies to a redshift that is larger than the high- 
redshift tail of our distribution, and we choose Zc = 0.525 
(made to coincide with one of the redshift nodes at which 
the k-|-E corrections are explicitly computed). 

There are many technical differen ces b etween the mod- 
els of Til and IMaraston et al.l l|2009l ) (e.g. IMaraston et all 
l2009l uses broadband optical colours whereas Til uses op- 
tical spectra; the fitting philosophies are different). For the 
current work, the primary difference is that Til uses the 
fossil record of galaxies at a given re dshift to infer the i r past 
star formation history. In contrast, IMaraston et all Hooi) 
fits a model to the colours of galaxies observed over a range 
of redshifts. To do so, they have to apply k+E corrections to 
all gala:xies in order to construct a sample that is potentially 
coeval. The risk of such an approach is that the inferences 
about the evolution of the sample reflect the assumptions 
used to define it, entered here via the k+E corrections - they 
only retained galaxies at a certain redshift that they predict 
would have been observed at any other redshift given the 
survey's selection function and the assumed stellar evolu- 
tion model, thus potentially considering a biased subsample 
of galaxies. In our approach, however, evolutionary colour 
paths that stray from the colour selection box are accept- 
able and crucially important for the results of our paper. 



In Til we assumed that that LRGs do not necessarily 
share the same stellar evolution or dust content, and they 
naturally allow for any intrinsic variation in LRGs, according 
to their luminosity, observed r — i colour and redshift. 

Brief ly, the Til mode l s are computed by running 
VESPA (|Toieiro et all |2007| . l2009l l on 124 high signal-to- 
noise stacks of LRG spectra, constructed based on the ob- 
served r — i colour, luminosity and redshift of the galaxies. 
VESPA returns a star-formation and metallicity history, as 
well as a present-day dust content, for each of the 124 stacks. 
These stellar evolutionary paths are turned into a spectral 
energy distribution, which can be computed at any redshift 
larger than the redshift of the stack and provide any set of 
colours, magnitudes, stellar masses or rest-frame luminosi- 
ties. 

An added subtlety lies in modelling the evolution the 
surface brightness of the galaxy, which is needed to fully 
characterise the survey selection function. This relies on 
modelling the physical size of a galaxy with redshift. Here 
we simply assume the physical size of the galaxy does not 
change with time, and we let the apparent size change with 
the angular diameter distance. In practice we know this is 
likely to be a simplification, as a growing body of literature 
discusses a potenti al size evolution of ea rly-type galaxies (see 
the Introduction of lSaracco et al.ll201ll and references within 
for a summary). These studies focus on the detection of a 
population of compact early- type galaxies at z; > 1, with 
effective radii a few times smaller than the effective radii 
of their local counterpa rts of the same stellar mass (e.g. 
Ivan Dokkum et "al]|2008l ). This does not seem to necessarily 
be a problem theoretically, as hydrodynamical simulations 
show that minor mergers can act to increase the effective 
radius of a galax y by a factor of a few between z = 3 and 
the present day (|Naab et al.ll2009l ). Nonetheless, the diffi- 
culty in characterising the (small) samples and making the 
measurements at high redshift, as well as the detection of a 
populatio n of compact early -type galaxies in the local Uni- 
verse (e.g. ICappellari et al. makes any realistic mod- 
elling of the size of LRGs with redshift unfeasible. Given the 
small redshift range of our sample and, crucially, the small 
number of objects that are discarded solely due to failing the 
surface brightness cut (< 0.01%), assuming no evolution in 
the physical size of our LRGs has a negligible effect on our 
results. 



4 MEASURED QUANTITIES 

The statistics that we use to test the evolution of the galaxies 
are the galaxy number a nd luminosity densi t ies, b oth as a 
function of redshift. As in lToieiro fc Percivall ( |2010l ). we use 
a proxy for the luminosity of the galaxy as L = 10~*^»/^'^, 
where Ali is K-l-e corrected to Zc = 0.525, and we work in 
arbitrary units of luminosity density throughout. 



4.1 The observed quantities 

We compute 



Vaz ^ Co 



(5) 
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sSAz 



(6) 



where Cg is an estimate of the spectroscopic completeness at 
the location of each galaxy, and corrects for the small frac- 
tion of galaxies in th e target sample that were not observed 
l|Percival et al.ll2007l ). Az is the width in redshift in which 
n{z) and £{z) are computed. We sample n{z) and £{z) in 
around 50 bins, between Zmin = 0.2 and Zmax ~ 0.45. 



4.2 The predicted quantities 

We begin by defining a redshift range, at low redshift, that 
we use to predict the number and luminosity densities to 
higher redshifts. This redshift bin, which we denote Azbasc, 
contains galaxies with redshifts in an interval can be varied 
in order to learn more about the evolution of LRGs. For 
each redshift greater than those in this bin we compute 

1 1 



■ E 



-Vobs,g{z) 



(7) 




0.20 



0.25 



0.30 0.35 
redshift 



0.40 



0.45 



■ E 



'-Vobs,giz). 



(8) 



4.3 Rates of change 

We define a number loss rate as 



nobs{z) 



(9) 



where At = t{Az) — t{Aztase), and a luminosity loss rate as 



riiz) 



ed(Z) 



(10) 



At V ^06^(2:) 
We can compute the average luminosity lost per galaxy 



(11) 



This is also the best estimate for the true merger rate of the 
galaxies that are the progenitors of the galaxies in Az^ase, 
assuming that any contaminants at redshift z have the same 
luminosity distribution as the progenitors. 



5 RESULTS 

In this Section we show our results for a selection of Az^ase 
and magnitude ranges of local galaxies. If the full sample 
was made of a coeval population of galaxies changing these 
values would have no effect on the derived merger rates. 
From previous work we know this assumption is likely to 
fail on two accounts: the targeting is less robust aX z < 0.2 
as can be seen by the fact galaxies at these redshifts tend to 
have small amounts of recent to interm ediate star formatio n 
and show an upturn in dust content (|Toieiro et al.l I2OIII ): 
and the fainter obj ects show stronger evid ence of a non- 
zero merger history (jToieiro fc Perciva]||2010l ). As our stellar 
evolution models are independent of the selection function 
to larger redshift, each of these two effects should become 
apparent as we vary the magnitude and the redshift of our 
local sample. Our interpretation is presented in Section [S] 



Figure 1. The observed and the predicted number and luminos- 
ity densities (first and second panel from the top, respectively). 
Number density in units of 10~* Mpc~^/i^ and the luminosity 
density in arbitrary units and Mpc~^?i^. The observed values are 
in black, and the predicted values are in red for the number and 
blue for the luminosity densities. The shaded area shows the red- 
shift range used to compute the predicted values. The third panel 
shows the loss rates in numbers (red) and luminosity (blue) given 
by Eqs. [9] and 1101 Bottom: the galaxy luminosity loss rate given 
by Eg. 111! All rates are in Gyr~^. 



We begin by making no selection in luminosity, and 
comparing observed and predicted densities from galaxies 
with Azbase = [0.20, 0.22] in the top two panels of Fig. [T] 
On the data side, we see a roughly constant number den- 
sity from z ~ 0.23 to 0.33, an excess relative to this plateau 
at 2: < 0.23, a bump aX z ~ 0.34, and a steep decline on- 
wards. The roughly constant densities are a direct result of 
the targeting algorithm, which aimed to select LRGs up to a 
fixed number density by following a passively evolving stel- 
lar population. The fact that these are roughly constant can 
in itself put some constraints on th e dynamical evolution of 
the sample (|Wake et al.lbOOdliooi '). The excess at low red- 
shift is likely to be contamination to the sample, where the 
target algorithm has less distinguishing power. The bump at 
z « 0.34 is a result of the widening colour cuts at that red- 
shift, as we move from cut I to cut II in the target selection. 
The steep decline after that is dominated by the apparent 
magnitude cuts, with the decline in the observed number 
and luminosity densities being steeper than that predicted 
from the low redshift data. 

The loss rates given by Eqs. [51 and 1101 are shown on the 
third panel of Fig. [T] and we combine the information in 
all four observables to compute the merger rate, or galaxy 
luminosity loss rate, given by Eq. 1111 which we show on the 
bottom panel of the same Figure. The merger rate is roughly 
constant and of the order of 3-4%, up to the point where the 
sample is dominated by the apparent magnitude limit. 

In Fig. O we show how the predicted quantities change 
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if we instead compute them from galaxies in different ranges 
Aztase- Note that differences in the predicted quantities (in 
the blue and red lines) across the different plots provides 
information about the nature and evolution of the galaxies 
at each of the redshift intervals - this is discussed further in 
Section [51 

We can also investigate how the loss rates change with 
the magnitude of the galaxies. As all absolute magnitudes 
are K+e corrected, a fiat cut in Mi with redshift selects 
the same population of galaxies in the absence of mergers. 
In Fig. [31 we show how the observed and predicted quan- 
tities change when we restrict the analysis to galaxies with 
Mi < —23: this corresponds roughly to the top 30% bright- 
est galaxies in the LRG sample. On the data side, we see 
many features in the redshift distributions disappear, lead- 
ing to a flatter curve. The exception is at low redshift, where 
the excess in density is still present. We show the same red- 
shift ranges of Aztase as explored in Figs. [T] and [2l Once 
again the observed curves are the same across all four pan- 
els, with differences in the red and blue curves holding infor- 
mation on the evolution of bright LRGs. We interpret and 
discuss these results in the next Section. 

5.1 Errors 

The errors quoted throughout this paper are Poisson, and 
assume that our ability to compare predicted with observed 
quantities depends only how many objects we have available 
to do so. We account for the change in photometric errors 
using Eq. [31 which is folded into the predicted quantities. 

It is particularly hard to estimate the effect of errors 
from the recovered star-formation histories on the computed 
merger rates. In TIf we presented an extensive analysis of 
the errors (see their Section 3.4), which we summarise here 
for completeness. There are two sources of errors that can 
affect the recovered star formation histories from VESPA: 
photon noise, which perturbs the spectra in way that we 
assume we understand; and systematic errors in models or 
parametrization, which are the result of our imperfect ability 
to correctly model certain stages of stellar evolution, chem- 
ical enrichment or dust extinction. The former is easily es- 
timated, via a Mon te Carlo approach that we explain in 
iToieiro et all (|200i). In this case we are in the regime of 
extremely high signal-to-noise spectra and photon error is 
negligible. In order to understand systematic errors, in Til 
we estimated these using the residuals of the best fit solution 
to the data (see their Fig. 3) - we found these errors to be 
at least 10 times larger. We use this error vector to obtain 
the SFHs used in the current analysis. 

The effect of the uncertainties in the SFH on the mea- 
sured rates comes from how many galaxies leave the survey's 
selection window as a function of redshift, and how bright 
they are. The window function of the survey is designed 
to follow a passive evolution, so only star-formation histo- 
ries that depart significantly from this model will affect the 
rates. An episode of recent to intermediate star formation of 
anything larger than 5% by mass is enough to take galax- 
ies out of the survey box for up to and around IGyr (the 
exact number depends on the length of episode, the Stel- 
lar Population Synthesis (SPS) models used, the metallicity 
and the exact age of the burst) and, even though none of our 
formal error estimates allow for a discrepancy of this size, 



solutions obtained with different SPS models can differ by 
a few times this amplitude - we explicitly address the effect 
of using different SPS models in our results in Section [6.41 

The effect on the merger rate in turn depends on how 
many galaxies are removed from the sample. Our assumption 
in Tll is that all galaxies within a cell follow the SFH given 
by the stack that represents that cell. Each of our Azbase 
samples is typically represented by 4 to 12 stacks (depend- 
ing on the luminosity range of the sample), and if one stack 
predicts a significant star- formation event (> 5% by mass), 
then 1/12 to 1/4 of the galaxies will be predicted to leave the 
sample for a significant period. This means that the effect 
on the merger rates are somewhat quantised, rather than a 
smooth function of an error in the SFHs. The effect can ob- 
viously be quite dramatic, and at its extreme it can predict 
that almost all galaxies would leave the sample, suggesting 
that the observed red sequence is an intricate combination of 
galaxies becoming red and blue throughout cosmic history 
(see Section [Ol for an example). 



6 DISCUSSION 

6.1 Interpreting the observables 

If the sample of SDSS-II LRGs was not affected by mergers, 
and assuming the SFHs are correct (see Section [5.1|) . then 
the observed luminosity and the number densities would ob- 
viously match the extrapolated values. Any deviations from 
this, implies merger activity and/or a non-zero net fiow of 
galaxies into the sample. We cannot determine a balance 
between mergers into and out of the sample. However, by 
comparing number and luminosity evolution we can start to 
understand the evolution that we do see. 

Fig. [4l shows schematically the simplest interpretation 
for three scenarios: 

(a) ripred < riobs and ipred = iobs - the simplest inter- 
pretation is that galaxies within the sample have merged 
together to decrease the predicted number density of ob- 
jects, whilst retaining the sample luminosity. This assumes 
no luminosity loss to the intra-cluster medium. 

(b) shows how contamination from objects at a given red- 
shift could also act to give ripred < nobs in the same way 
as internal mergers. By contaminants here we refer to ob- 
jects that would enter the sample for a short time only, and 
do not evolve to be present in the sample at low redshift. 
However, the comparison between lobs and Ipred helps to 
differentiate the two scenarios, as in this case we expect the 
excess in the observed number density due to contamination 
to be matched by an excess of luminosity density. The two 
excesses will be matched in amplitude provided the contami- 
nants have the same luminosity function as the galaxies that 
evolve through the sample. 

(c) Upred > nobs and Ipred > lobs- TMs scenario requires 
an enrichment of the sample towards lower redshift - ei- 
ther by the merging of fainter galaxies becoming brighter 
and crossing over the faint magnitude cuts, bluer galaxies 
turiung redder and crossing over the colour cuts, or a com- 
bination of both. 

Note that in the case of fainter galaxies merging together to 
become brighter and entering the sample by crossing over 
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0.20 0.25 0.50 0.55 0.40 0.45 0.20 0.25 0.50 0.55 0.40 0.45 0.20 0.25 0.50 0.55 0.40 0.45 



Figure 2. The same as Fig.[T] but exploring different intervals for Azfj^se^ used to compute the predicted quantities. 
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predicted 
high-z 



observed 
low-z 



observed 
high-z 



a) 



-O^ 



internal mergers: Tlobs > f^predj ^obs — ^pred 



b) 



■O^ 



contamination: riobs > T^pred-,^obs > ^pred 







c) 

external mergers: Tlohs < T^predi (-obs < (-pred 



Figure 4. Diagram of possible interpretations for discrepancies in 
tile observed and predicted quantities. Tlie middle column shows 
the observed luminosity and number densities. The solid arrows 
show the prediction given by the stellar evolution from VESPA 
and the survey's window function, and they lead to the left hand 
side column that shows the predicted luminosity and number den- 
sities. Note that whereas VESPA can only predict changes in the 
luminosity density, the window function can also affect the num- 
ber densities. The window function is fully known and introduces 
no ambiguity when computing the predicted quantities - in here 
we show the numbers being conserved just for simplicity. The 
right hand side column shows the observed number and lumi- 
nosity densities, and the dotted arrows represent our proposed 
interpretation. See Section l6.ll for more details. 



the magnitude cuts, scenarios (a) and (c) represent the same 
physical process - the only difference between is where the 
survey's magnitude cut lies. We explore this fact in a little 
more detail in Section [6.31 

Obviously, we can expect multiple scenarios such as 
those in Fig. |4]to happen simultaneously. If this is the case, 
under the assumption that the luminosity function of the 
contaminants matches that of the objects in the sample, and 
that luminosity is conserved in a merger event. Eg. Illl gives 
the best estimate for the merger rate required for galaxies 
within the sample between and Azbase- In fact, even 
if these assumptions are broken, Eq. [11] still acts as a lower 
limit on the true merger rate. This is mostly easily seen from 
Equation (|lip - we are effectively under-estimating Ipred, as 
some luminosity has been lost to the intracluster medium, 
resulting in an under-estimation of rg. 



6.2 Interpreting the results 

Using the results of Section O we can identify redshift inter- 
vals that are likely to be dominated by each of the scenarios 
shown in Fig. jj] Starting with the results for all LRGS, 



shown in Fig. [TJ we see the predicted number and luminos- 
ity densities both fall short of the observed values up to 
z « 0.38, suggesting a combination of scenarios (a) and (b). 
The bump at z « 0.34 is not reproduced by the stellar evo- 
lution of present-day galaxies, strongly suggesting that is it 
caused by contamination of galaxies as the colour cuts widen 

- in other words, the wider breadth of colour in the sample 
at those redshifts is not predicted by the stellar evolution of 
more local galaxies. Therefore the signal here is likely dom- 
inated by scenario (b) . Finally, the sample of local galaxies 
predicts a much shallower slope in the number densities as 
we go towards higher redshifts, suggesting that they were 
still being assembled (in a dynamical sense) at those epochs 

- in other words, had the local galaxies still be in one piece 
at 2 > 0.4 they would have been bright enough to be in the 
sample. This puts us in a regime likely to be dominated by 
scenario (c) at these high redshifts. 

Moreover, at the high-redshift tail where scenario (c) 
dominates, we note that the number loss rate is higher (in 
amplitude) than the luminosity loss rate. In other words, 
the typical luminosity per object at 2 > 0.38 is lower than 
expected. This suggests that the LRG growth due to ex- 
ternal merging is happening predominantly at the fainter 
end of the sample. We will investigate this hypothesis fur- 
ther in Section [6.2.21 Note that if this growth was constant 
with luminosity then both loss rates would have a negative 
value, but would be matched in amplitude - resulting in a 
zero galaxy luminosity loss rate, interpreted in this case as 
a zero internal merger rate. 



6.2.1 Changing Azftase 

We now consider the effect of changing the low redshift inter- 
vals Aztasc, from which we infer the galaxy evolution from 
high redshift as shown in Fig. (2] Changes in the inferred 
loss and merger rates as we change our low-redshift sample 
are sensitive to differences in the low redshift population of 
galaxies. 

We find no significant difference in the evolution of 
galaxies observed between z — 0.23 and z = 0.27 - see the 
first two panels of Fig. [2] - but comparing the rate of lumi- 
nosity loss per galaxy (lower panels in the figures) to those 
in Fig. [T] suggests that the lower redshift galaxies experi- 
ence a more pronounced growth. This is consistent with a 
hypothesis in which the low redshift end of the sample is 
not as pure as the rest of the sampl e - the targeting algo- 
rithm is less robust at these redshifts l|Eisenstein et al.]|200ll : 
iToi eiro et al.ll201ll ). Furthermore, the sample will naturally 
be more enriched towards fainter galaxies, which could have 
a different growth rate - see Section [6.2.21 for more details. 

Similarly, galaxies observed at redshifts between 0.33 
and 0.35 are seen to have a different growth compared with 
galaxies at lower redshifts. This is not unexpected, as the 
bump in number and luminosity density at « 0.34 is not 
predicted by galaxies at low redshifts. These results suggest 
that the wider colour selection of cut II selects galaxies not 
only with a different stellar evolution, but also a different 
merger growth rate - essentially a different population of 
galaxies. For the same reasons as before, we argue that this 
growth must be happening predominantly at the faint end. 
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6.2.2 Changing the luminosity 



In Fig. Owe restrict our analysis to galaxies with a K+e cor- 
rected Mi < —23, which corresponds roughly to the bright- 
est third of the galaxies in the sample. We show results for 
the same four ranges of Azbase as before. It is immediately 
clear that a lot of the features in the observed quantities 
are less pronounced - the bump at z ~ 0.34 almost dis- 
appears, and the high-redshift tail flattens out as we re- 
strict ourselves to what is close to magnitude-limited sam- 
ple (colour cuts excluded). The exception is the bump at 
low redshifts, which, combined with the measured loss and 
merger rates for Azbase = [0.20, 0.22], suggests the luminous 
galaxies observed at low redshift are still growing, and this 
growth is approximately matches that of lower luminosity 
objects. Clearly this also fits with the hypothesis that the 
lower redshift galaxies form a population with more dynam- 
ical evolution than the higher redshift galaxies within the 
sample. 

For all other cases of Azbase we find a picture that is 
consistent with passive dynamical growth, with the only - 
and marginal - evidence for external growth coming for the 
galaxies in Azbase = [0.33,0.35], perhaps suggesting that 
not even the brightest galaxies are purely coeval across the 
two LRG colour cuts described in Section [S] 



6.3 Probing the typical luminosity of mergers 

We noted in Section fG.ll that scenarios (a) and (c) represent 
the same physical merger process, but can be distinguished 
as (c) caused objects to cross the apparent magnitude cut on 
the sample. In Section [6.21 we suggested than (c) dominates 
at high redshift, while (a) dominates at low redshift. We can 
use the information available to define a typical luminosity 
for the dominant merger process happening. 

For each case of Azbase we take the redshift at which 
the number and luminosity loss rates cross the zero line and 
compute the absolute magnitude that would correspond to 
the apparent petrosian magnitude equal to the survey's limit 
of Tp — 19.5 - in other words, the faintest absolute magnitude 
observed at that redshift. These are the galaxies that are 
beginning to enter the sample when dynamical growth tran- 
sitions from external to internal. Note that the two growth 
rates change sign at very similar epochs, so we can use ei- 
ther without significant effect on our results. In Fig. [5] we 
show how this transitional absolute magnitude changes for 
the four different cases of Azbase (those shown in Figs.[l]and 
[2J. 

Once again our result is that the dynamical growth is 
happening predominantly at the faint end, for galaxies ob- 
served in all the four intervals of Azbase we studied. 



6.4 Robustness to SPS models 

In lToieiro et~al] l|201lh we showed how the measured stellar 
evolution of LRGs depends on the SPS models used - see Fig. 
6 in that paper for a summary of the differences. The results 
presented in this paper thus far assume the stellar and colour 
evolution given by t he Flexible Stellar P o pulation Synthe- 
sis (F SPS) models of IConrov et"al] (|2009D : ]Conrov fc GunnI 
which gave better fits to the spectral data and re- 
sulted in the more passive stellar evolution model in our 
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Figure 5. Tlie dots sliow the transitional absolute magnitude as 
a function of Azbase (ou the right-hand side y-axis) - see Sec- 
tion 16.31 This is the absolute r — band magnitude that corre- 
s ponds to the = 1 9.5 cut on petrosian magnitude of cut II 
l|Eisenstein et alJlioOlh at the redshifts for which we infer a tran- 
sition from external to external growth. For context we also show 
the distribution of absolute magnitudes of galaxies in the sample 
(on the left-hand side y-axis). 



previous study. However, as we could not conclusively de- 
cide which set of SPS models gives a more accurate descrip- 
tion of the real stellar evolution of LRGs, in this section 
we discuss how the stellar evolution given by other models 
would change the interpretation given above. The results on 
the recovered merger rate and interpretation of the dynam- 
ical evolution can be dramatically affected by events of star 
formation. 

In the case of the Maraston et al. 2011 models (MIO, 
submitted), the amount of recent to intermediate star for- 
mation dominates only in galaxies at 2 < 0.3. The resulting 
interpretation is that contamination of the LRG sample ex- 
tends to a larger redshift than what would be assumed using 
the FSPS models, and that these local galaxies are not the 
evolutionary products of the LRGs observed at higher red- 
shifts. At 2: > 0.3, however, the recovered dynamical evolu- 
tion is qualitatively similar to that obtained with the FSPS 
models, although with slightly larger rates. 

The interpretation using the stellar evolution given by 
the models of Bruzual & Chariot (200:^) (BC03), however, 
is incompatible with the hypothesis that the LRGs form a 
coeval population of galaxies since z < 0.45. The significant 
and continuous (with redshift) amount of recent to interme- 
diate star formation lead to a consistently under-prediction 
of the number and luminosity densities at high redshifts from 
more local samples, in a way that is only explained by a con- 
tinuous amount of contamination in the sample. 



6.5 Comparison with previous results 

The resu lts in this paper ar e in qu alitative agreement with 
those in iToieiro fc Percivall (|2010l '). where we used the lu- 
minosity weighted power-spectrum to argue that the LRG 
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population is growing via external mergers at the faint end. 
Quantitatively, however, the rates we find in this paper are 
lower by roughly a factor of two. The main difference be- 
tween the two appro aches lies on the treatment of the stellar 
evolution - whilst in Toicir o fc Perciv al (20f0) we assumed 
that t he passive stellar evolution model of Maraston et al] 
(|2009l ) was suitable for all galaxies in the sa mple, in here 
we ins tead use the stellar evolution models of iToieiro et al.l 
which use the fossil record of LRGs of different colour, 
luminosity and redsh ift to infer the most - likely stellar evolu- 
tion. Furthermore, in iToieiro fc Percivall ||2010D . we assumed 
that the scatter in colour and magnitude was due to intrin- 
sic or photometric errors, which is a byproduct of assuming 
that all LRGs in the sample share the same stellar evolution. 
Instead, the analysis presented in this paper makes no such 
assumption, and is robust to an influx and/or an outflux of 
objects if we interpret the measured rates as a lower limit. 
It is the independent treatment of the stellar evolution, and 
the detailed comparison of predicted with observed quanti- 
ties that allows us to make an assessment on the nature of 
the mergers without the need for a clustering analysis. 

Comparisons with other works is less trivial, due to the 
typically different samples used. We refer the reader to Sec- 
tion 8 of Tojeiro fc Percival ( 2010. ) for a summary of recent 
results and comparison. We continue to argue that merg- 
ing at the faint end of the LRG population is a more likely 



explanation published for the published evidence for both 
luminositv growth ( Brown et al. 20071: IWhite et al.l l2007l: 
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7 SUMMARY AND CONCLUSIONS 

In this paper we have demonstrated how we can use the fossil 
record of galaxies to extrapolate from low redshift samples 
to match samples selected at higher redshift. By compar- 
ing against the observed number and luminosities, we can 
measure growth and consider where the dynamical evolu- 
tion of galaxies is required. We have applied this technique 
to consider the dynamical evolution of galaxies in the SDSS- 
II LRG sample, although the technique could be used to link 
different surveys at low and high redshift. 

For the SDSS-II LRG sa mple, we have p reviou sly pre- 
dicted the stellar evolution in IToieiro et al.l (|201ll ). by us- 
ing the fossil record encoded in the spectra of LRGs. This 
allowed us to compute the most likely stellar evolution of 
observed LRGs in a way that was decoupled from the sur- 
vey's selection function. We used this information to make 
predictions on the number and luminosity of LRGs at higher 
redshift from a local sample of galaxies, and we showed how 
differences in these quantities can be interpreted as a mini- 
mum merger rate. Our main conclusions can be summarised 
as the two following points: 

• the LRG population is not purely coeval, with some of 
galaxies targeted at 2; < 0.22 and at z > 0.34 showing dif- 
ferent dynamical growth than galaxies targeted throughout 
the sample; 

• the LRG population is still dynamically growing, and 
this growth must be limited to the faint end; 



Based on these results, we argue that a coeval popu- 
lation of galaxies that has been dynamically passive to less 
than 1-2% can be robustly selected from the LRG sample, by 
imposing a cut on the K-|-e corrected absolute magnitude, 
and by selecting galaxies with z G [0.23,0.45]. This merger 
rate is lower than those previously measured by roughly a 
factor of two, because, by using empirically determined evo- 
lutionary tracks we no longer need to assume passive evolu- 
tion, and this previously led to incorrect evolutionary tracks 
for some galaxies. 

We caution against interpreting these results beyond 
the limits of the stellar population models, and once again 
stress the important of good models and statistically good 
fits for analyses of this type. However, given the power in- 
herent in the type of approach, and improvements in SPS 
modelling due in the future, we consider this analysis to be 
a forerunner of significant future work. 
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